Effect of impurities on supersolid condensate: a Ginzburg-Landau approach 
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Basing our arguments on a wave function that contains both positional and superfluid order, we 
propose a Ginzburg-Landau functional for a supersolid with the two order parameters necessary to 
describe such a phase: density nB(r) and supersolid order parameter ipv- We argue that adding 
lighter ^He atoms to a ^He supersolid produces attractive regions for vacancies, leading to patches of 
higher Tc. On the other hand, the supersolid stiffness decreases in this granular state with increased 
^He disorder. Both effects are linear in "^He concentration. 



INTRODUCTION 



WAVE FUNCTION FOR A SUPERSOLID 



Recent experiments by Chan and Kim Q have gen- 
erated renewed interest in the possibility of the existence 
of a supersolid phase in ^He. There is pioneering theo- 
retical work by, among others, Andreev and Lifshitz Q, 
Reatto 0, Chester j^, Leggett |^ and Anderson 0. Re- 
cently, Anderson |^ and Anderson, Brinkman and Huse 
1^ have discussed the problem. Recent developments 
have been reported at a KITP workshop and are available 
online 

The main purpose of purpose of this note is to de- 
velop an approach that allows a discusssion of the effect 
of ^He impurities on a ^He supersolid. In what follows, 
we show how the approach leads to a Ginzburg-Landau 
functional for the supersolid. We propose to follow the 
collective degrees of freedom that are candidates for the 
supersolid properties, namely the vacancies that are in- 
trinsic to an incommensurate solid. Thus, we follow 
Andreev and Lifshits |^, and more recently Anderson 
1^, who argued that the supersolid requires some defects, 
which are inherently present in the phase at T = 0. Va- 
cancies and interstitials are obvious candidate defects in 
the ^He solid. For simplicity, we focus here on vacancies 
In this picture, the lattice sites are equivalent, there 
is no disorder and yet the He solid is not commensurate. 

Thus, we can consider the ground state of ^He at T = 
as a state in which there is a non-zero concentration of 
vacancies. Although a localized vacancy may require a 
non-zero activation energy, the fact that they are mobile 
(the more so due to the large zero-point motion of the 
He atoms) can bring the bottom of their band to zero 
energy. For any non-zero concentration ny of vacancies, 
they will form a Bose condensate at T = 0. However, to 
properly describe the quantum nature of the condensate 
we have to address not the densities but the quantum 
amplitudes of the relevant degreees of freedom. Thus 
the Ginzburg-Landau (GL) functional for the supersolid 
deals with the vacancy Bose field i/jy = ^/ny ex.p{ia) as 
an order parameter. 



In this section, we discuss wave functions that display 
solid order with, at the same time, a Bose condensed 
fraction of vacancies. 

We start with a lattice state of hard-core bosons; the 
boson occupancy cannot exceed one per lattice site. The 
operator creates a ^He atom at site i of the lattice. 
The commensurate crystal wave function is then 



|M/o) = ]^5+|0) = | T,T,...T), 



(1) 



corresponding to a "ferromagnetic" state, with no vacan- 
cies. For some purposes, it is convenient to use a spin-1/2 
representation 12]: 



bi, St = b+h - 1/2. 



(2) 



To introduce vacancies, a simple procedure would be to 
apply Ny times to create Ny vacancies. However, 

it has been understood since BCSjind was pointed out in 
the present context by Anderson 8] that a phase coherent 
superfluid state must be a linear combination of states 
of different particle number. To achieve this, we apply a 
tilt operator to the "ferromagnetic" commensurate state. 
The rotation operator in the spin representation is (now 
using Pauli operators) 

]^exp(io-i • ki(l)i/2) = Y[[cos ]-(t)i + i{(Ti ■ ki) sini0i].(3) 



Here, (j)i is the tilt angle at the z-th site and is the 
rotation axis, which lies in the x, y plane. The latter's 
angle with respect to the y axis will be denoted by a^. 
Since ± z. 



(Ti ki = {Si a. +5^ a+). 
With af = iay = ±zexp(±zai), we get 



(4) 



(5) 
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We operate with Eq. (3) on |^o) = Ili^z^lO) ^^e 
Eq. (5) to find the tilted state 

I*) = l[[cos^<Pi + i{ie'"*)sm^<p,bi]bt\0) 

i 

= l[[cos^cj>,bt-e"^*sm^cj,m, (6) 

i 

where we have used the fact that because of the hard-core 
constraint, = hf gives zero acting on |^o)- The state 
1^) is a hnear combination of states with specified phases 
and different numbers of vacancies; it is normahzed and 
the expectation value of vacancy occupation is: 

nV = 1 - (5+6,) = sin2 (7) 

The vacancies are mobile due to rearrangements of the 
atoms, which is facilitated by their large zero-point mo- 
tion. We model the kinetic energy of vacancies with a 
simple nearest-neighbor hopping hamiltonian Its 
expectation value in |^) is |l4j : 

KE = -t^ibfb.^h.c.) 

(ij) 

= (— 1/2) ^sin(/)^ sin(/)j cos(a^ — aj). (8) 

(ij) 

With all ai the same, |^) is a state with both solid 
order and a phase coherent vacancy contribution. The 
vacancy order parameter is 

V'F = (^.) = ^sin(/.e-'". (9) 

In the spin language, when all are the same, |^) 

is a state with maximal total spin S = M/2, where M 
is the number of lattice sites. The expectation value of 
= Y.S! is {S') = (M/2) COS ^. 

GINZBURG-LANDAU FUNCTIONAL 

We start with with the postulate that the GL func- 
tional of a supersolid phase in ^He has to contain two or- 
der parameters: density ^^(r) and superfluid amplitude 
V^y(r). Although there are only ^He atoms, these atoms 
particpate in two distinct phenomena: solid and super- 
fluid. In the region of the experimental phase diagram 
0, where one enters a possible supersolid phase from 
the solid at T ~ lOOmK at flxed pressure, solid order 
is presumably already well established and hence is well 
outside of the GL regime. Thus, we concern ourselves 
with a possible second-order transition across a normal 
solid to supersolid transition line. The amplitude [^/^^(r)! 
is small at the transition and is in the GL regime. The 
wave function |^) in Eq. © offers a unified description 



of both solid order with periodically modulated density 
nB(r) and vacancies with the order parameter i/jy plus 
the constraint n^(r)+ny(r) = 1. In the discussion below 
we focus on the tpv field. 

Toward developing the GL approach, we assume that 
the vacancy density is small (as is consistent with the 
experimental situation) so that the tilt angle (f) is small. 
To estimate the superfiuid stiffness, we examine the KE 
as the phase a acquires a slow variation from site to site. 
From Eqs. (8,10)and small 0, we get 

AKE={t/ao) J ^^r|V^yp(Va)^ (10) 

where ao is the lattice spacing in a simple cubic lattice. 
We see that the superfiuid stiffness is determined by the 
hopping amplitude t and the density of vacancies 

Next we examine the effect of a potential that couples 
to the ^He density. Hi = Uih^hi. From Eq. (7), 

^ J dru{r)\^v{r)\', (11) 

where Q is the unit cell volume, u{r) = Q UiS{r — Vi) 
is the potential energy and we ignored the constant 
J dru{r). In the last line, we took advantage of the 
approximation that (j) is small so that, from Eq. (10), 
^ ny = sin^ ^(j). In our proposed state |^), a po- 
tential that is repulsive for ^He atoms would be attractive 
for vacancies, tending to increase their local density. 
Our discussion leads to the total GL free energy density 

m = [a(;^-l)-^i(r)/^2]|^y(r)p + (^/ao)|^y(r)|2|Vap 

^ C 

+ b\i^v{v)\\ (12) 

The quartic term in F represents the restriction against 
putting two vacancies in the same place, as well as dy- 
namical interaction terms. 

EFFECT OF DISORDER 

Chan et al |2] find the remarkable result that by adding 
^He atoms in the ^He solid, the measured Tc grows and 
the superfiuid stiffness ps drops as a function of ^He con- 
centration. The data is reproduced here in Fig. 1. 

Let us therefore address the effect of disorder on the su- 
persolid within the GL framework. In particular, because 
of its experimental relevance, we are interested in the ef- 
fect of ^He impurities on the superfiuid properties. One 
might expect that disorder would suppress the supersolid 
state as it would interfere with the phase ordering as it 



3 



E 




XJppm] 



XJppm] 



FIG. 1: Dependence of transition temperature and superfluid 
stiffness on "^He concentration. Data of Cfian et. al. |3, also 
quoted in ^ 

does in superconductors. However if tfie supersolid arises 
due to tfie presence of vacancies, tfien tfie presence of de- 
fects tfiat create vacancies or tfiat can increase tfie local 
concentration of vacancies could lead to enhanced super- 
fluidity. Here, we offer some speculations on the effects of 
a non-zero concentration of ^He atoms on the transition 
temperature and superfluid stiffness. 

1) The presence of lighter ^He atom defects increases 
the zero point fluctuations locally; this repels nearby host 
atoms. As we have seen, a repulsion for "^He atoms is an 
attraction for vacancies; this increases the local vacancy 
density as seen in Fig. 2. This leads to an increase in 
Tc. To first approximation, the effect would be linear in 
the ^He concentration. 2) At the same time, there would 
be an effect on the superfluid stiffness p^. We suggest 
that ps would decrease with increased ^He concentra- 
tion because ^He sites would produce attractive sites for 
the superfluid condensate and localize it more. While 
^He creates puddles of enhanced supersolid condensate, 
the stiffness as a reflection of the global phase rigidity 
would drop on average with increased ^He doping. In 
other words, to the extent that substitutional impurities 
inhibit the exchanges of ^He atoms that are required to 
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enable vacancy hopping, we expect that any foreign atom 
substitution will diminish the effective leading to a re- 
duction in superfluid stiffness. Again, the effect will be 
linear in impurity concentration, see Fig(^). 

So we propose that a ^He site is attractive for vacan- 
cies and we model the effect of ^He substitution as a 
random attractive potential for vacancy density with the 
potential u{T)/ft = -1-Ap3(r), where A (> 0) is the cou- 
pling strength of the effective ^He - ^He repulsion. The 
corresponding contribution to the GL free energy density 
is 



-An3|^y(r)|^ 



(13) 



Here, we have taken an annealed average, (p3(r)) = 
as the average density of ^He atoms. 

We see, from Eqs. (12-14), that the effect of ^He sub- 
stitution is to increase the transition temperature from 
TO to 

re = (l + ^n3)T°, 5TJT!=^ (14) 

an effect that is linear in ^He concentration ns. One could 
estimate the parameters by examining the experimental 
data. In fact the numbers from Chan et al l2j show an 
enormous effect. From the data. Fig. [3 we find that 
shift 6Tc = 150 mK, or STc/Tc ^0.75 for an increase in 
^He concentration Sns of about 10 ppm. Put differently, 
the as defined critical temperature nearly doubled upon 
adding 0.001% impurities. To date, no known superfluid 
states exhibit such a stong sensitivity to inipurities. This 
is one of the many puzzles in this system |l5j . 

Now we turn to stiffness corrections. We argued above 
that the stiffness term will soften in the presence of ^He 
defects. The effective hopping t in the GL free energy 
density of Eq. (13) is then t = to{l—gns/nv) with ^ > 0. 
Since ps cx we conclude that 

Ps=p^si^-gns/nv), p^, = 2{to/ao) \^v\^ . (15) 

From the data, we estimate Sps/p^ ^ 1/4 for Sns ^ 10 
ppm. Assuming ny ^ 0.1%, we find g ^ 25. 



CONCLUSION 



Vacancy attractive potential near 3He site 
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FIG. 2: When included in the ^He lattice, ^He atoms will 
push the host atoms away. This redistribution is equivalent to 
an attractive potential for vacancies, thus generating a larger 
vacancy density near ^He. 



Although we motivated our arguments with a wave 
function appropriate for considerations in the grand 
canonical ensemble, in our discussion of the effect of ^He 
impurities on the superfluid stiffness, we assumed implic- 
itly that changes in the average concentration of vacan- 
cies nv with ^He concentration are much less important 
than the decrease in the kinetic energy. We have no ev- 
idence for this but strictly speaking we cannot exclude 
the possibility that adding ^He increases the equilibrium 
vacancy density. If this is the case, we need to augment 
our arguments by explicitly allowing ny be dependent on 
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723 . This would be an extension of our analysis that could 
lead to modifications of our results. Such an increase in 
nv would be a simple source of the Tc increase with ad- 
dition of ^He. However, even the presence of patches of 
higher Tc could increase the bulk Tc by proximity effects. 
Similarly, the mass decoupling that is measured in the 
torsional oscillator would have higher onset temperature 
due to patches of higher Tc. We leave this issue until the 
experimental situation is clearer. 

We presented a GL approach to the supersolid phase 
that contains two order parameters, density modulation 
n5(r) and superfluid ?/;y(r). Both the solid and super- 
fluid fields are of course made from the same ^He atoms 
and realize a dual quantum mechanical behavior of these 
atoms. As others have argued, it seems reasonable to us 
that lattice defects, for example vacancies, are crucial for 
supersolid phenomena. 

We find that disorder from ^He atoms produces regions 
of depleted ^He density, thus attracting vacancies. This 
attractive potential due to disorder can be though of as 
creating locally regions of higher Tc in the GL functional 
for t/jy. At the same time the granularity in ipy produced 
by disorder leads to suppression of superfluid stiffness. 
Both of these effects are linear in concentration of ^He 
atoms. 
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